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Sequential suppression of quarkonia and high-energy nucleus-nucleus collisions 
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According to the usual application of the sequential-suppression picture to the dynamics of heavy 
quarkonia in the hot medium formed in ultrarelativistic nuclear collisions, quark-antiquark pairs 
created in a given bound or unbound state remain in that same state as the medium evolves. We 
argue that this scenario implicitly assumes an adiabatic evolution of the quarkonia, and we show 
that the validity of the adiabaticity assumption is questionable. 

PACS numbers: 25.75.Nq, 12.38.Mh, 14.40.Pq 



More than 25 years ago, finite-temperature gauge- 
theory studies on the lattice of the screening of static 
(color) charges prompted Matsui and Satz to suggest that 
in a quark-gluon plasma (QGP) the formation of bound 
charmonia may be prevented [1]. It was quickly realized 
that the various cc or bb states might in fact be disso- 
ciated at different temperatures 0,13]. Over the years, 
paralleling progress in theoretical studies of quarkonia 
properties (see Refs. 0, [I| for recent reviews), this led 
to the "sequential-suppression" picture of heavy quarko- 
nia as QGP thermometer Q. According to the latter, a 
given state will be totally suppressed above a threshold 
temperature — which might actually be smaller than the 
transition temperature to a QGP. 

This prediction is supported by several approaches. 
First, spectral functions are extracted from lattice-QCD 
computations of correlators of quantum numbers for 
heavy quark-antiquark pairs [7H10|. The disappearance 
of a peak in the spectral function then signals the sup- 
pression of a state. Alternatively, one resorts to an effec- 
tive in-medium quark-antiquark (QQ) potential — cither 
derived from lattice-QCD computations of spatial corre- 
lators (U-H, fldl-fljl or derived within finite-temperature 
field theory }16l4l8l | — , which then enters a Schrodinger 
or Bethe-Salpeter equation, whose bound states model 
the quarkonia in the medium. The suppression of a given 
state takes place when it is no longer bound by the poten- 
tial, although the precise criterion for dissociation might 
be open to discussion (lij . 

In either description, the in-medium bound states of 
heavy quark-antiquark pairs are, be it explicitly stated 
or not, eigenstates of a Hamiltonian. Note that as the 
effective potential might actually possess an imaginary 



part [15l - ll7l |. these eigenstates are not necessarily sta- 



ble, but might have only a finite lifetime. That notwith- 
standing, the generally accepted picture is that of a 
temperature-dependent suppression pattern, in which at 
a given energy density of the medium some states sur- 
vive, while more excited ones are not bound and thus do 
not form. 

This picture is seemingly supported by experimental 
measurements in ultrarelativistic nucleus-nucleus colli- 
sions. At the SPS, the NA50 collaboration reported that 



the anomalous suppression of the ip' in Pb-Pb collisions 
sets in at a smaller average in-medium path crossed by 
the cc pair than for the J/ip [1(3]. At the much higher 
LHC energy, the CMS collaboration studied bottomonia 
and observed yields consistent with the idea that the ex- 
cited T(2S) and T(35) states are more suppressed than 
the deeper bound T(15) 0. 

The usual explanation for such observations in high- 
energy nucleus-nucleus collisions is the following, where 
for the sake of simplicity we leave aside so-called "initial- 
state effects" At an early stage after the collision, say 
some instant to, the created deconfined medium reaches 
high enough an energy density that a given quarkonium 
state, which we shall refer to as "excited" , is suppressed, 
while another state of the same system, hereafter the 
"ground state", is bound. The common lore is then 
that, as the medium expands and cools down (t > to), 
the ground state stays unaffected, whereas the depopu- 
lated excited state remains suppressed, even when the 
medium temperature has dropped below its dissociation 
threshold. The only possibility left to the excited state 
for being recreated is at the transition to the hadronic 
phase, through the "recombination" of till then uncorre- 
cted heavy quarks and antiquarks 22 , 23| . That scenario 
constitutes the standard implementation of the idea of 
sequential suppression of heavy quarkonia in heavy-ion 
collisions. Note that this description totally ignores the 
possible finite lifetime of the ground state due to the 
imaginary part of the in-medium potential. 

Inspecting the scenario sketched above critically, it re- 
lies on two basic ingredients. There is first the sequential- 
suppression pattern in the "initial condition" at to, whose 
theoretical foundations we discussed above. The second 
element in the scenario, which to our knowledge has not 
been examined before, is the implicit assumption that 
"the quarkonium ground state remains the ground state" 
over the duration of the medium evolution. Recasting 
this statement more mathematically, a QQ pair initially 



1 When comparing relative yields of different states of a given sys- 
tem, say S-channel charmonia or bottomonia, for a fixed type of 
nuclear collisions, these effects should play a minor role. 
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in the eigenstate with lowest energy of the (effective) 
Hamiltonian describing in-medium quarkonia remains in 
the lowest-energy eigenstate. More generally, the same 
will hold for every initially bound state — up to late elec- 
troweak decays which take place outside the medium. 
That is, it is assumed that heavy quarkonia are contin- 
uously evolving eigenstates of an adiabatically changing 
instantaneous Hamiltonian. Accordingly, the scenario for 
the sequential suppression of quarkonia in the medium 
created in high-energy nucleus-nucleus collisions relies 
on the hypothesis that the effective in-medium quark- 
antiquark potential varies slowly enough that each QQ 
pair is at every successive instant in an energy eigen- 
state. We now wish to investigate the validity of this 
assumption. 

Before going any further and to dispel any confusion, 
let us note that the adiabaticity we discuss in this note 
is neither that of the medium evolution — related to the 
production of entropy — , nor the adiabatic assumption a 
la Born-Oppenheimer which allows one to separate glu- 
ons from the nonrelativistic heavy quarks when writing 
down an effective potential for the latter 24] . 

Let \n(t)) denote the eigenstates of an instantaneous 
Hamiltonian H(i), with respective energies E n (t). Fol- 
lowing the approach of Ref. [,2S] , a common criterion for 
the validity of the adiabatic theorem is the requirement 
that for every pair of states \n(t)), \n'(t)) and at every 
instant in the evolutional 



amplitude for a matrix element of dV/dT between eigen- 
states of the instantaneous Hamiltonian 



\(n'(t)\H(t)\n(t))\ 
[E n (t)-E n ,(t)] 2 
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with H(i) the time derivative of the Hamiltonian. In the 
case of interest in this note, H(t) coincides with the time 
derivative dV/di of the effective QQ potential. In turn, 
the latter is simply the product of the rate of change T 
of the medium temperature times the derivative dV/dT 
of the potential with respect to T, where for the sake of 
simplicity we have assumed that the medium is (locally) 
thermalized. For T, we took the results from a simulation 
of central Pb-Pb collisions at the LHC within dissipativc 
hydrodynamics 26] , considering the evolution of temper- 
ature at the center of the fireball: within the first 7 fm/c 
of the evolution (that is, as long as T > 200 MeV), f 
always remains larger than about 30 MeV per fm/c and 
up to 50 MeV per fm/c in the early stages. 

For the QQ potential, we considered the lattice QCD 
results of Ref. [27j , including the parameterization 



V(r) 



-A^/l+N f /6Tg 2loop (T)r 



(2) 



with A ~ 1.4 and g2\oop(T) the 2-loop perturbative cou- 
pling, where 0.5 < a s (T) < 1.0. One then finds as typical 



, , dV 

n (*) n{t) 



200- 500 MeV- fin. 



We use a system of units in which h = c = 1. 



The numerator in Eq. (JTJ is thus of order (80-160 MeV) 2 . 
In turn, the denominator is of order (100-350 MeV) 2 for 
the excited bb states, so that the ratio can be in some 
cases smaller than 0.1, for other channels larger than 1. 
Because of those channels, it is far from warranted that 
the adiabaticity assumption holds: the potential evolves 
so quickly that a quark-antiquark pair which at some 
time is in a given instantaneous eigenstate will a short 
while later no longer be in the evolved eigenstate, but will 
have components over all the new eigenstates — including 
the new ground state, which shows that even if crite- 
rion ([T]) holds for the latter, yet it is populated by con- 
tributions from excited states. 

We wish to emphasize here that this "repopulation" 
mechanism is neither the customary recombination at 
hadronization, nor the feed-down from late decays, but 
a natural consequence of the "reshuffling" of QQ states 
due to the rapid medium evolution. 

A naive picture of the effect of this rapid evolution is 
provided by dividing the typical size r rms w 0.3-0.75 fm 
of a bound bottomonium by the characteristic velocity 
v ~ 0.3c of the nonrelativistic constituent quark and an- 
tiquark, which gives a duration r rs 1-2.5 fm/c for an 
"orbit" of the b quark. On such a time scale, the QGP 
cools down by 30 to 75 MeV, resulting in a significant 
change in the effective potential ([2]), which illustrates 
why the adiabatic evolution of bottomonia is far from 
being warranted. 

As a final argument against using the hypothesis of an 
adiabatic evolution of QQ pairs in a QGP, we note that 
recent studies emphasized the fact that even when cri- 
terion ([1]) is satisfied — i.e., the evolution is slow — , the 
system with evolving Hamiltonian can be driven from 
one instantaneous eigenstate to a different one at later 
times by resonant interactions (see e.g. Ref. [28]). The 
latter lead to Rabi oscillations between eigenstates — that 
is, they are tailored to induce transitions which violate 
the adiabatic theorem — on a time scale given by the in- 
verse of the Rabi frequency ujr. 

In the case of a QQ pair in a quark-gluon plasma at 
the temperatures found in high-energy nuclear collisions, 
there are obviously plenty of degrees of freedom around 
with energies matching possible transition lines. The cor- 
responding Rabi frequencies however depend on the in- 
teraction term. Adopting for the sake of illustration a 
dipolar interaction, one finds values of ir/ujR, which in a 
two-level system is the time after which a transition has 
occurred with probability 1, of the order 2 to 20 fm/c, 
depending on the transition Bohr frequency, the medium 
size and the assumed coupling strength. This means that 
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on such a time scale a QQ pair certainly does not re- 
main in the same instantaneous eigenstate, which again 
hints at the invalidity of the adiabatic theorem for heavy 
quarkonia in a dynamical QGP. 

One might be tempted to argue that in an effective- 
potential approach, the transition-inducing degrees of 
freedom have been integrated out. Yet the construction 
of an effective theory ultimately relies on the adiabatic 
theorem (29[, so that it is inconsistent to use the notion 
blindly here. More precisely, we surmise, although we 
have not investigated this idea in detail, that the viola- 
tion of adiabaticity caused by resonant interactions trans- 
lates into the imaginary part of the effective in-medium 
potential, which physically has the same effect of giving 
a finite lifetime to the Hamiltonian eigenstates. 

In summary, we have shown that the usual scenario for 
the sequential suppression of heavy quarkonia in the hot 
medium created in ultrarelativistic heavy-ion collisions 
relies on the hypothesis that QQ pairs evolve adiabat- 
ically in the medium. We have then presented several 
arguments which make us doubt that this assumption 
holds. In our view, this hints at the idea that a given 
quark-antiquark pair is in a constantly changing linear 
combination of instantaneous energy eigenstates, rather 
than in a smoothly evolving unique eigenstate. As a con- 
sequence, one should explicitly follow the evolution of the 
pair in the QGP, using a dynamical microscopic descrip- 
tion, as attempted in Refs. [I(J|3l[. Such an approach 
should of course be suited for rapidly evolving media, so 
as to eventually be able to compare with experimental 
results from high-energy nucleus-nucleus collisions. On 
the other hand, in the regime of a (quasi-)static environ- 
ment of the quarkonia it should also make contact with 
equilibrium-based formalisms as lattice QCD or finite- 
temperature field theory. 
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